Universal Linear Density of States for Tunneling into the Two-Dimensional Electron 

Gas in a Magnetic field 
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A new technique permits high fidelity measurement of the tunneling density of states (TDOS) of 
the two-dimensional electron gas. The obtained TDOS contains no distortions arising from low 2D 
in-plane conductivity and includes the contribution from localized tunneling sites. In a perpendicular 
magnetic field, a pseudogap develops in the TDOS at the Fermi level. Improved sensitivity enables 
resolution of a linear dependence of the TDOS on energy near the Fermi energy. The slopes of this 
linear gap are strongly field dependent. The data are suggestive of a new model of the gap at low 
energies. 
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Characteristics of electrons tunneling into and out of 
a two dimensional (2D) system differ considerably from 
those of the three dimensional (3D) case. The distinction 
is especially pronounced when a magnetic field is applied 
in the tunneling direction (perpendicular to the plane in 
the 2D case). In the simplest picture, such a magnetic 
field effectively localizes electrons in the 2D system. An 
electron tunneling into an energetically unfavorable posi- 
tion cannot readily move away and instead tends to move 
in circles. As a result, tunneling measurements of 2D sys- 
tems in a magnetic field display effects attributable to a 
"pseudogap" in the tunneling density of states (TDOS) 
at low injection energies [B-k§. In contrast, for a 3D sys- 
tem the tunneling electron can move parallel to the field 
lines to evade being localized at a position of high po- 
tential energy and such a field-induced gap has not been 
detected. 

There have been two chief obstacles in interpreting 
tunneling measurements of a 2D electron system. First, 
the in-plane conductance must be kept much larger than 
the tunneling conductance, otherwise the measurement 
represents transport within the 2D plane instead of across 
the tunnel barrier. This is a major problem when the 
2D electron density is low or when the 2D electron sys- 
tem acts effectively as an insulator in the quantum Hall 
regime. Second, schemes such as tunneling between two 
2D layers [g]||] yield a convolution of effects from both 
layers. Tunneling from a 2D into a spectroscopically fea- 
tureless 3D layer achieves superior resolution of 2D fea- 
tures 0. However, until now it was only possible to 
measure zero-bias tunneling between 2D and 3D in semi- 
conductor heterostructures M. 

In this letter, we report results from a new technique 
which we call "time-domain capacitance spectroscopy" 
(TDCS) for measuring the I-V characteristics of struc- 
tures to which direct ohmic contact is not possible. Us- 
ing TDCS, we measure the tunneling current into a 2D 



electron system at arbitrarily low values of in-plane con- 
ductivity. Unlike other tunneling current measurements, 
our technique detects all of the current entering or exit- 
ing the 2D layer, including those arising from electrons 
entering localized sites. 

Zero bias suppressions in tunneling conductance are 
well known to be a many body phenomenon. Differ- 
ent approaches have been taken to solve this complex 
problem involving electron-electron interactions, disorder 
and magnetic field [pHPl. Our result indicates that this 
field-induced tunneling suppression differs qualitatively 
from previous theoretical predictions. First, the TDOS 
is found to have a universal linear dependence on energy 
near the Fermi energy for all field strengths and electron 
densities. The slopes of this linear gap are strongly field 
dependent. Second, the high excitation tunneling spec- 
trum shows a change in curvature as the field strength is 
increased. 

Figure la depicts the type of samples used in our ex- 
periment. The 2D electron system is sandwiched between 
two electrodes, close enough only to the bottom electrode 
to permit tunneling of electrons. Mesas etched from two 
wafers grown using molecular beam epitaxy have been 
studied. Both wafers consist of a degenerately n doped 
GaAs substrate followed by an AlGaAs tunnel barrier. 
On top of that a GaAs quantum well is grown which de- 
fines the 2D electron system. A thick AlGaAs blocking 
barrier prevents charge transfer between the well and the 
top GaAs electrode. The blocking barrier contains an n 
doped region to provide electrons for the well. The first 
wafer (wafer A) 1 1 has been studied in detail previously 
using frequency dependent capacitance measurements JM 
to determine zero-bias tunneling characteristics. The sec- 
ond wafer (wafer B) |L2| has a smaller dopant concentra- 
tion and a thinner tunnel barrier. A DC bias applied 
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FIG. 1. (a) Structure of our MBE grown samples, (b) Evo- 
lution of the conduction band profile of our sample in one 
measurement cycle, (c) External circuit including a standard 
capacitor and a HEMT used to measure the current flowing 
out of the sample. 

to the top gate permits variation of the density of elec- 
trons in the quantum well from depletion to 6 x 10 n cm -2 
(sample A) and 3 x 10 11 cm -2 (sample B) . The mobility 
of the 2D electron system in sample A was estimated O] 
to be 1 x 10 5 cm 2 /Vs at a density of 2 x 10 n cm~ 2 , and 
the mobility of sample B is expected to be considerably 
higher. TDCS enables us to measure away from zero bias 
and extract the complete I-V characteristics. All features 
of the tunneling suppression described in this paper have 
been observed in both samples. 

Figure lb shows the evolution of the conduction band 
diagram for our samples during one cycle of the TDCS 
measurement. We start with the 2D electron system in 
equilibrium with the 3D substrate. At time t = 0, a 
sharp (< 10ns step rise) voltage step is applied. This cre- 
ates an offset in the chemical potential on the two sides 
of the tunnel barrier, inducing a tunneling current. As 
electrons tunnel, this offset equilibrates, and the decay 
signal is recorded in real time. To measure the current 
across the tunnel barrier, a capacitance bridge is utilized 
(Fig. lc). Voltage steps of opposite polarity are applied 
to the substrate of the sample (Cr) and to one plate of a 
standard capacitor (C$). The other plate of C$ and the 
top electrode of Ct are electrically connected, and the 
signal at this balance point is fed into the gate of a high 
electron mobility transistor (HEMT). Through a 70MJ1 
resistor Rbias the DC bias of the HEMT is established. 
The amplitude of the step applied to Cs is adjusted so 
that the voltage at the balance point decays to zero after 
electrons cease tunneling. All measurements described in 
this letter take place at a temperature of 65 mK. 

Before the voltage step is applied, the 2D plane has 



the same electrochemical potential at all positions. Im- 
mediately after the voltage step is applied (t = + ), no 
charge has been transferred into the 2D electron system. 
The voltage pulse is not screened by the quantum well 
and the 2D plane remains an equipotential. The simple 
planar geometry of the sample dictates that this remains 
true even in situations of very low 2D conductivity. At 
t = + , the voltage across the tunnel barrier (Vbarrier) 
is simply a fraction of the voltage step applied (V s t e p), 



given by V b 
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with a and b defined in 



Fig. lb. The 3D substrate is always highly conducting, 
and electrons can tunnel everywhere into the 2D electron 
system. 

Using charge conservation we have determined that 
there is a fixed relationship between the barrier current 
and the displacement (measured) current. Immediately 
after the applied voltage step, the current across the bar- 
rier is found to be proportional to the initial time deriva- 
tive of the voltage at the balance point: 



dV b 



^barrier 



= C^ where C s = c^-c.Cs-CCs (1) 



The capacitances C\ and Ci in equation (1) are the sim- 
ple geometrical values C\ = (A/(b — a) and C2 = eA/a, 
where e is the dielectric constant and A is the area of the 
mesa. Elsewhere ["JJ, we prove that the relationship (1) 
remains true independent of thermodynamic DOS vari- 
ations in the quantum well. By applying voltage steps 
of different amplitude and taking initial time derivatives 
of the corresponding transistor signal, the complete I-V 
characteristics of the tunnel barrier can be mapped out. 
Signals from the experiment are extremely faint, and im- 
mense signal averaging is involved in our measurements. 
Each point on an I-V trace may require averaging of as 
many as 100,000 time traces. A novel signal processing 
and rapid averaging system |lfj] permits data acquisition 
with enormous (~18 bit) digital resolution. 

Figure 2 shows the tunneling conductance (I/V) of 
sample A plotted against the voltage across the barrier 
for magnetic field strengths of 0, 1, 2, 8 and 16 Tesla at a 
fixed electron density of 1.9 x 10 n cm -2 . This density is 
high enough so that no zero-bias tunneling suppression 
is observed at zero field. Application of a magnetic field 
reduces the tunneling conductance around zero bias. The 
suppression becomes deeper and wider as the field is in- 
creased. This field-induced tunneling suppression differs 
quantitatively from the logarithmic suppression [|16[ in 
the low density regime at zero field [O. Near zero bias, 
the conductance is found to have a universal linear de- 
pendence on the excitation voltage for all magnetic field 
strengths and electron densities. Moreover, an increase 
in the strength of the suppression is accompanied by a 
change in the curvature of the high excitation part of the 
conductance curves when the magnetic field is increased, 
as shown by the bottom inset of Fig. 2. Even though the 
conductance curves at high excitation appear 
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Voltage across tunnel barrier(mV) 
FIG. 2. Dependence of the tunneling conductance (I/V) 
of sample A on excitation voltage across the tunnel barrier 
for different magnetic field strengths at 65 mK and a fixed 
density of 1.9 x 10 11 cm~ 2 . Bottom inset: same set of con- 
ductance curves, zooming in near zero bias. Top left inset: 
Conductance curve for sample B at 4T and v ~ 3, showing 
Landau levels at higher excitations in addition to the zero 
bias suppression. Top right inset: Simulated I-V for tunnel- 
ing between two 2DEG's, using 16T data from sample A (see 
text). 

rounded at high fields, the zero bias region remarkably 
remains linear in voltage, with both the magnitude and 
the slope significantly reduced. This singular behavior is 
illustrated by the insets of Fig. 3, which zoom in near the 
zero bias region of the conductance curves for different 
field strengths. Such a linear energy dependence of the 
TDOS is observed over the full range of densities in both 
samples except near depletion (n < 5 x 10 10 cm -2 ). The 
top left inset of Fig. 2 displays a conductance curve from 
sample B at a field of 4T and v ~ 3. In addition to the 
zero bias suppression, features associated with adjacent 
Landau levels can be identified at higher excitations. 

In order to compare our data to results from double 
well experiments |2jJ^] , we compute the I-V curves ex- 
pected for tunneling between two 2D electron systems: 

reV 

J oc / g{E- eV)g{E)dE 
Jo 

Here both 2D systems are assumed to have identical 
TDOS g(E) deduced from our 2D-3D tunneling data 
from sample A at 16 Tcsla. The resulting I-V curve, 
as shown in the top right inset of Fig. 2, qualitatively 
resembles that from double well experiments. 

To our knowledge, no existing model other than the 2D 
Coulomb gap p^% predicts such a universal linear DOS 
at low energies for such a wide range of field strengths. 
However, contrary to expectations for a simple Coulomb 
gap the slopes of the observed linear gap are strongly 
field dependent. Figure 3 shows the slopes of the linear 
regions of the conductance curves plotted against inverse 
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FIG. 3. Insets: Tunneling conductance of sample A as a 
function of excitation voltage for 6 different magnetic field 
strengths. The units are the same in all insets, with the 
abscissae in mV and the ordinates in /i-Siemens. Different 
ranges are chosen to display the singular behavior near zero 
bias. Main figure: Dependence of the slope of this linear gap 
on inverse field strength. 



magnetic field strength. For filling factors v < 1, the 
data points fall nicely onto a straight line extrapolating 
to a negative intercept on the vertical axis. For low fields, 
there are deviations from the straight line as the filling 
fraction varies between integer and non-integer values. 

In the Coulomb gap picture, the states in the vicinity 
of the Fermi level are assumed to be localized. These 
electrons are treated as classical point charges at fixed 
positions in space with no overlap of the electronic wave- 
functions. The phase space available for electron tun- 
neling is reduced since it costs more energy to add an 
electron to the system when another electron is located 
close to the tunneling electron. The resulting Coulomb- 
gap in the TDOS varies linearly with excitation energy in 
2D with a slope determined solely by physical constants 
such as the electronic charge and dielectric constant |L7j . 
In contrast the slope of the linear gap in our data depends 
inversely on field strength with an offset, suggesting that 
a simple Coulomb gap is inadequate in explaining the 
tunneling suppression in our experiment. 

We propose a phenomenological model which describes 
the data quite well. This model is inspired by one pre- 
viously developed for tunneling into a system of random 
sized metal particles |18J. In this picture, the 2D sys- 
tem is modeled as isolated puddles with uniform charg- 
ing energies and random background offsets. Interac- 
tions among the puddles are neglected, in contrast to the 
Coulomb gap model. This assumption may be justified 
due to the presence of the nearby 3D conducting sub- 
strate which screens the interactions among the puddles. 




FIG. 4. (a) From the top: I-V curve and corresponding 
conductance curve of a single puddle. A random offset volt- 
age shifts conductance curves of different puddles along the 
voltage axis. The bottom figure shows the V-shaped con- 
ductance curve resulting from summations of these randomly 
shifted curves, (b) Puddles of different sizes produce conduc- 
tance curves with different slopes, giving rise to an overall 
conductance curve shown at the bottom, (c) When the pud- 
dles become small enough that e/C becomes larger than the 
background offset voltage, the overall conductance will vanish 
near zero bias, (d) In this regime, when puddles of different 
sizes are taken into account, a U-shaped conductance curve 
will be obtained. 

Even with no applied voltage, the Fermi level of each 
puddle will not, in general, align with the Fermi level of 
the 3D bottom electrode. The energy of a puddle with 
capacitance C increases or decreases by e 2 /C when an 
electron is added to or removed from the puddle. Thus 
the Fermi energies on the two sides of the tunnel bar- 
rier are only equilibrated to within e 2 /2C of each other. 
The Fermi energies in these puddles are assumed to be 
uniformly distributed within this range because of, for 
instance, a random background voltage offset created by 
nearby dopants or impurities. Each puddle contributes 
a Coulomb blockade type FV characteristic, leading to 
a conductance curve which is constant in voltage except 
for a region of width e/C randomly displaced from zero 
bias where the conductance vanishes (Fig. 4a). 

The sum of conductances from all puddles will thus be 
linear in voltage near zero bias. The slope of this linear 
gap is inversely proportional to the capacitance C, while 
the width of the gap is proportional to C. Since the ca- 
pacitance C is proportional to the area of the puddles, 
our data can be explained if the average area of the pud- 
dles varies inversely with magnetic field strength. This 
model assumes that the high voltage conductance of an 
individual puddle is proportional to its area and the total 



area occupied by puddles is constant. 

Another appealing feature of this "Coulomb blockade 
gap" model is that it can explain the different curva- 
tures of the tunneling spectrum at high excitations as 
well as the negative offset of the slope vs. inverse field 
strength dependence shown in Fig. 3. In the low field 
limit, the range of random energy offsets is larger than 
the Coulomb blockade energy. When puddles of differ- 
ent sizes are taken into account, the resultant TDOS will 
be a superposition of linear gaps with different widths 
and slopes, giving rise to a negative (convex) curvature 
at high excitations while preserving the linear behavior 
at zero bias (Fig. 4b). 

In the high field limit, some of the puddles become 
small enough so that their Coulomb blockade energies 
exceed the range of the background offset energy. Their 
conductance curve will no longer be V-shaped. In this 
regime, conductance contributions from puddles of a par- 
ticular size will be zero at low bias up to a certain voltage 
beyond which the overall conductance rises linearly with 
voltage to the unsuppressed value, as depicted in Fig. 4c. 
Summing contributions for puddles of various sizes leads 
to the U-shaped conductance curve in Fig. 4d, concave 
with respect to voltage at high biases. As long as there 
exist some puddles large enough to produce a V-shaped 
conductance curve, the linear behavior of the overall con- 
ductance is preserved near zero bias, albeit with a much 
reduced slope. This argument can be carried further to 
explain the finite magnetic field required to produce a 
zero slope in the TDOS as extrapolated from our TDCS 
data. This happens when the puddles are small enough 
so that the Coulomb blockade energy of every puddle in 
the system exceeds the range of the background offset en- 
ergies. It is not necessary to have an infinitesimal puddle 
to achieve a zero slope for conductance near zero bias. 

The above deductions are based on the assumption 
that larger puddles break up into smaller ones and that 
the mean area of the puddles shrinks linearly with in- 
creasing field strength. Electrons therefore charge par- 
allel plate capacitors whose lateral dimension is propor- 
tional to the magnetic length. From the width of the gap 
in our data, we estimate the proportionality constant to 
be about 6. We do not yet have a clear answer to the 
question of what the puddles are or why they shrink as 
the magnetic field strength is increased. 

While this simple picture of a Coulomb blockade gap 
may not provide a complete description of the system, it 
seems to be able to explain qualitatively most features 
in our data. A more thorough understanding of the tun- 
neling suppression will require inclusion of interaction, 
charging effects and the quantum mechanical properties 
of the 2D electron system. 
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